NASA-CR-204731 

I^T MYS ' CAI J ° URNAL ' 469:L 61-L64, 19% September 20 

' T1 ’ e Amc,K;an Astronomical Society. All rights reserved. Primed in U S A. 




& 

helioseismic constraints on the gradient of angular velocity at the 
base of the solar convection zone tLUUTY AT THE 

w u/ Urv r- A. G. Kosovichev 

■ ■ ~ *— — 

ABSTRACT 

rotation of the convTcUon^ne “‘"Th interi . or ‘° the latltud,nal differential 

sphttmg coefficients of thep-mode frequencies, obtained ^durinc lSfi 3 ^ ° R fth " Sun ' Usi "M rotational 

we have found that the thickness of the transitional laver is n no + n na P° at ! be Bl ^ Bear ^ olar Observatory, 
of the transition occurs beneath the ad, aba Sly sSfied oao ” d ".<‘ 3 * 28 “"»• a " d ttat most 

dynamo theories of the 22 yr solar activity cycle * part of the convection zone, as suggested by the 

^ Ct headmSS ' Sun-rotatIo7 meth ° dS '’ ^ ana ' ySiS ~ Sun: aCtivity ~ 1 Sun: interior ~ Sun: oscillations - 


1 INTRODUCTION 

Helioseismology has established the existence of a layer of 
■ rong gradients of angular velocity at the base of the solar 
convection zone (e.g., Brown et al. .989; Goode et al 1991 
Tomczyk, Schou, & Thompson 1995). This layer separates the 

Z'f rb 0nC r CXhibiting stron 8 latitudinal differLtial rota- 

n from the radiative interior rotating almost rigidly. Turbu- 
lence generated in the layer is likely to mix material in the 
upper radiative zone, resulting in the observed deficit of Li and 
Be (see, e.g., Zahn 1992). However, the theoretical estimates 

C, ,rL» r r C ' S ”> lf U,0n | and thC dneknessofthe transition layer 
( tachocl.ne ) depend on details of turbulent energy and 
momentum transport, and are uncertain (Spiegel & Zahn 

Perhaps of the greatest interest, the transition layer is the 

w°th- P aCC f ° r the S ° ]ar dynamo <scc - e -g- Weiss 1994) 
Within this layer, the toroidal magnetic flux that appears at the 

surface m various forms of solar activity is generated from the 
radial component of the poloidal field (Branderburg 1994) 
The toroidal flux is believed to be mainly accumulated in a thin 
layer just beneath the convection zone because convection 
would quickly destroy the toroidal flux if the layer were widely 
extended into the convection zone. However, as recently 

3thffi y h 8er ,b & Brandcnh urg(1995), this layer cannot be 
() ^' h " b ^ tause thc P enod of the solar cycle, which depends 
llrhulc ' nt magnetic diffusion time through the laver 

0 05R ~ ~ G/ hCy C u‘ imated thc thickness “» be at least’ 

th ^ i i ‘ where R is the solar radius and H is 

the local pressure scale height. p 

tr ; ,n^r at f ° f , the thickness and precise location of the 
ansition layer by standard helioseismic inversion techniques 
rom rotational splitting of oscillation p-modes are rather 
uncertain. Attempts to resolve the layer under global smooth- 
ness constraints lead to either an oversmoothed angular 
velocity proffle or to spurious oscillations around the transition 

st?ated C th?t°e de Ct f al ;J 991) - Thom P son (1990) has demon- 
forS h a • ,f hC trans,tion were discontinuous the 

Jeaffin Tn aSCISmiC mV p rSIOnS Sti " P roduce 3 broad smooth 
region To overcome these difficulties, Goode et al (1991) 

assumed that the transition layer is, in fact, discontinuous and 


found that their model fits the helioseismic data best when the 

In Z j"? COm . c,des with the base of the convection zone. 

IPhr rnf ’ We dem ° ns 5 ate that the discontinuous model of 

wffh a fn°r IS Ti Tl. m t0 the data and that a model 
with a transitional layer of finite thickness (*0.09 ± 0 04R) 

fits the data more accurately than the discontinuous model 
he midpoint of this layer is found at 0.692 ± 0.005R, slightly 
below thc convection zone. b y 

2 . FREQUENCY SPLITTING AND SOLAR ROTATION LAW 

Observed p-mode rotational frequency splitting is tradition- 
1986) CPreS m ^ (DuVa "’ Harvey ’ * Pomerantz 


= L J %H(n- /)/>*,, m 

* * 0 . 1 . 2 .... \LI 


(i; 


where n, /, and m are the radial order, the angular degree, and the 
‘ agalar order of a normal mode, respectively; L = \l(I + l )l> -• 

aTe seN a of fficSd T hus ’ the observational data 

L odd "-coefficients,” , , , for mode multiplets 
/ * F ° r thc modes of intermediate and high degree / w hich 

^dt e ,h con Tr !om - ,he ®“ dSJ 

reorns n l ^ funC f t,0n of the solar station law. 

order l PGnl -°£ aS ^ 0 iatcd ^"dre functions of 

’ k (^) (Kosovichev 1988), and expressed as 


Hfr, 6 )/ 2 tt = ^ ak A 


*-0.1.2.. 


» 


sin 6 ' 


( 2 ) 


T™n^^rlV ,2 ‘ ,(2k + '>"• The "**■ *«-« 


I K 

a k (n, l ) = 


AM 


X Wlj + /(/ + \)Vh - 2U„ , V„ / - V; u ] p H {3) 

where f/„, ( and V , are the radial and horizontal displacement 
eigenfunctions of oscillation modes, p(r) is the density, and 
nJ \ Unl + w + l)Vn.i]pr^ dr is the mode inertia. The 
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contribution of the last two terms in the buckets of equation 
(5) resulting from the Coriolis effect, does n 
S^n/tticrefoic can be treated as a small perturbation. 

Then we introduce coefficients a as 


i t i t * ■ T 


where 


a' k (n , /) 


/) = a’ k (n, l) (1 Ci,/)> 


* A h (r)[Ul, + /(/ + IWhW dr • 


"n. / 


(4) 


(5) 


New coefficients aUn, /) represent the “seismic averages of 
the radial function A k , with the weighting function propor- 
ttoaUtohe mode en'ergy dens,,,. and C ^ — - 
describes the effect of the Coriolis force (Ledoux 195 U- 
For high-frequency p-modes, a’ k (n, l ) are essentia y un - 
tions of only one asymptotic variable: the angu ar p ase spe 
T /L or equivalently, the radius of the lower turning point r 
of 'ie modes (Gough 1984). This asymptotic behavior is 
demonstrated in Figure 1, in which filled circles show t e 

Sic averages to, a mode, of so, a, rotation 
Figure 2. We use this asymptotic property of the sphtt g 
coefficient for the analyses of the observational data by 
averaging (n, 0 over short intervals of r, . 

3. ESTIMATES OF WIDTH AND LOCATION OF 
THE TRANSITION LAYER 

We have used the splitting ^-coefficients (a, , a 3 , and la 5 ) ) of 
approximately 800 multiplets of / - 
ard & Libbrecht (1993) from 1986 and 
tions made at the Big Bear Solar Observatory (BBSO) The 
a -coefficients corrected for the Coriolis effect, a k (n, ), 
each year have been grouped and 

in 16 equal intervals of r„ between 0.5 R and 09 R each bemg 
0 025ft wide. Then the weighted averages for the y 
been computed. The final averages, («».))* U ’ 

are shown by the diamonds with the error bars 'n F.gure T The 
variation at the base of the convection zone is Particularly 
strong for a\. Therefore, we have used this coefficient t 
estimate the thickness and the location of the transition layer 
The absence at the base of the convection zone of a shar P 
variation of a’, , which represents the seismic average of the net 
angular momentum of a spherical shell, means essentially zero 
net torque 1 * between the radiative and convection zone 
(Gough^l985; Gilman, Morrow, & DeLuca 1989). We sha 
discuss the limits of the variation of a) in a 

Helioseismic inversions have revealed that the radiative 
interior essentially rotates rigidly. Therefore .n romt.on law 
(2), functions A k (r ) for k > 1 ,aie approximately •**£*££ 
in the radiative zone. Also, the averaged aria, l) coeflraente 
shown in Figure 1 suggest that A, is almost constant n he 
convection zone. Therefore, we have considered A } in 
parametric form, 


A y (r) = A x () <t>(r), 


( 6 ) 
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N 
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where 4>(r) = 0.5(1 + erf [2(r - r 9 )/w]), erf is the error 
ffinction ro is the radius of the central point of the transition 
zone and w is the characteristic thickness corresponding to the 
variatkm of 4>(r) from 0.08 to 0.92. The value of A,„ 


440 



h. 



1 to, £» to .1 <«. 0 ■»— 

the BBSO 1986-1990 data (/ - 5-60). 

(*22 nHz) is determined from the flat part of ® » JJV ** 
convection zone. Then, we computed (a^(n,l )>mod coefficient 

to eu,X”(6) using «q»*i°" < 5 > and the Mm f 
procedure that was applied to the BBSO data, an , na y, 
determined the mean square difference, 

X 2 (W, r 0 ) = 2d T lX a 3j)ol» _ (U 3j)roodl • (7) 

j=i 

The function *>, r.) f shown in Figure 3, reaches ^a mini- 
mum of 16.8 at r, * 0.692 R and w * 0.09 R. The boundaries 
of the shaded area correspond to the increase of x Y > 

Itr in r 0 and w. A similar result was obtained using different 
approximations for the transition function, 4 »(r) and by 
choosing different intervals of r, for averagmgthe da a Th 
analysis repeated with an extended BBSO data set 
?/ = 5-140) provided to us by P. R. Goode gave the same result. 
( We note ffiat the position of the center of the transition 
layer is determined more accurately from the data than is the 
thickness. It is evident that the center of the layer is beneath 
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Fig. 2. — Contours of I2(r, B)/2tt from 328 to 461 nHz, separated by 7 nHz 
of a solar rotation model used in Fig. 1. Dashed curve shows the lower 
boundary of the convection zone. In this model, based on the first three terms 
of eq. (2), A, (r) is approximated below the convection zone by a constant and 
it is approximated in the convection zone by two linear functions, the first of 
which was chosen to fit gradually increasing a \ , while the second function 
(describing the subsurface shear layer) was found by matching the internal 
rotation with the surface rotation (Snodgrass 1992) — assuming the conserva- 
tion of the angular momentum in the subsurface layer (Foukal & Jokipii 1975). 
Estimated thickness of the subsurface layer is approximately 12 Mm (0.017/?). 
Functions Ay (r) and A${r) are taken in the form of eq. (ft). 


the boundary of adiabatic stratification of the convection zone, 
which is at 0.713 ± 0.003/?, and coincides with the zone of 
the sharp variation of the sound speed (Fig. 4). For Goode & 
Dziembowski’s (1991) model with a discontinuity at the base 
of the convection zone {w = 0, r„ = 0.71/?), \ I 2 « 36 is 
significantly higher than the minima! value. Therefore, their 



r 0 /R 

Fig. 3. — Contours of * 2 (>v, r 0 ) evaluated from eq. (7) at \ 2 - 18, 19, 21, 24, 
28, 3ft, 48, 64, and 9ft; w is the thickness of the transition layer, r 0 is its central 
radius. Shaded area corresponds to the increase of \ 2 by 1 from its minimum 
value, or 1 a uncertainty in the parameters. Error bars show the l <r uncertainty 
estimated from statistical modeling by adding Gaussian noise to the data. 


A 3 (nHz) r5u/u 



r/R 


Fig. 4.— Solid curve with the shadow (indicating 1 a uncertainty) shows the 
parameter Ay of rotation law (2), estimated from the BBSO data. Points with 
error bars represent the variations of the ratio of the pressure to the density, 
u = p/p, relative to a standard solar model, inferred from the GONG data 
(Gough et al. 1996). Vertical hatched column shows the location of the base of 
the adiabatically stratified part of the convection zone, as determined by 
Christensen-Dalsgaard, Gough, & Thompson (1991) and Kosovichev & 
Fedorova (1991). 


model with a very thin transition layer at the base of the 
convection zone can be excluded. 

4. CONCLUSION 

Estimated from the 1986-1990 BBSO solar oscillation data, 
the thickness (0.09 ± 0.04/?) and the central position 
(0.692 ± 0.005/?) of the transition layer from the nearly 
uniform rotation in the radiative interior to the differential 
rotation in the convection zone at the base of the solar 
convection zone are generally consistent with the require- 
ments of the aw-dynamo theory of the 22 yr solar activity cycle 
(Rudiger & Brandenburg 1995). Our results show that most of 
the transition occurs beneath the zone of adiabatic convection. 
However, the upper boundary of the transition layer extends 
into the convection zone. To determine whether the transition 
layer coincides with the zone of convective overshoot is an 
important problem of helioseismology. 

Another steep angular velocity gradient is likely to occur 
just beneath the solar surface because the angular velocity in 
the convection zone inferred form helioseismology is substan- 
tially higher than the angular velocity observed at the solar 
surface. Resolving the structure of the subsurface gradient 
layer is another challenging task for the new helioseismic 
projects, GONG and SOHO, and would require accurate 
observations of rotational splitting of high-degree p-mode 
frequencies. 
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